Trapping in scale-free networks with hierarchical organization of modularity 
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A wide variety of real-life networks share two remarkable generic topological properties: scale-free 
behavior and modular organization, and it is natural and important to study how these two features 
affect the dynamical processes taking place on such networks. In this paper, we investigate a simple 
stochastic process — trapping problem, a random walk with a perfect trap fixed at a given location, 
performed on a family of hierarchical networks that exhibit simultaneously striking scale-free and 
modular structure. We focus on a particular case with the immobile trap positioned at the hub node 
having the largest degree. Using a method based on generating functions, we determine explicitly 
the mean first-passage time (MFPT) for the trapping problem, which is the mean of the node-to- 
trap first-passage time over the entire network. The exact expression for the MFPT is calculated 
through the recurrence relations derived from the special construction of the hierarchical networks. 
The obtained rigorous formula corroborated by extensive direct numerical calculations exhibits that 
the MFPT grows algebraically with the network order. Concretely, the MFPT increases as a power- 
law function of the number of nodes with the exponent much less than 1. We demonstrate that the 
hierarchical networks under consideration have more efficient structure for transport by diffusion in 
contrast with other analytically soluble media including some previously studied scale- free networks. 
We argue that the scale-free and modular topologies are responsible for the high efficiency of the 
trapping process on the hierarchical networks. 
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I. INTRODUCTION 



Complex networks are a powerful and versatile math- 
ematical tool for representing and modeling structure of 
complex systems j and their wide applications in dif- 
ferent areas have made them become a subject of a large 
volume of research in the past decade [3, Within 
the general framework of complex networks, scientists 
can offer in qualitative terms the detailed microscopic 
description of structural properties and complexity of 
real-life systems. Extensive empirical analysis on diverse 
real systems has unveiled that many, perhaps most, real- 
world networks are simultaneously characterized by the 
two most remarkable features: scale-free behavior [5| and 
modular organization 6, 7, 8]. The scale-free nature of a 
network means that its degree distribution P{k) follows 
a power law as P{k) ^ k~'^ with the degree distribution 
exponent in the range of 2 < 7 < 3, while the modular or- 
ganization implies that the network is formed by groups 
(modules) of nodes that have a significantly higher in- 
terconnection density compared to the overall density of 
the whole network. The important finding of these two 
fundamental natures has led to the rising of research on 
some outstanding issues in the field of complex networks 
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such as exploringthe generation mechanisms for scale- 
free behavior [ll detecting and characterizing mod- 
ular structure [9l. [lOl. [Til , [l^ . and so on. On the other 
hand, it has been shown that the two characteristics are 
closely related to other structural properties such as av- 
erage path length [Hill and clustering coefficient [8|. 

In principle, one of the main reasons for studying struc- 
tural properties of complex networks is to understand 
how the dynamical processes are influenced by the un- 
derlying topological structure [lB|- Among a plethora 
of random processes, random walks with wide range of 
distinct applications to many science branches, have at- 
tracted a considerable amount of recent attention within 
the physics community [ip, [13, [H, [H, [13, [HI, [H, HI, [H, 
[H, [li, [13, [H, [M [M [3J Particularly, trapping issue, 
an integral major theme of random walks, is relevant to 
a variety of contexts, including target research |32j, i33|, 
photon-harvesting processes in photosynthetic cells [34| . 
and characterizing similarities between the elements of a 
database [35| has led to an increasing number of the- 
oretical and practical investigations over the last sev- 
eral decades. Numerous authors have made concerted 
efforts to study trappin g pr oblem in different media, in- 
cluding regular lattices '36| , Sierpinski fractals [13, , 
T fractal 13&|, small- world networks [40], and scale-free 
networks l4ll . [4p . [43l . [3, [1^, as well as other struc- 
tures [13, |47| . |48| |. These studies unclosed many unusual 
and exotic phenomena of trapping on diverse graphs. 
However, the trapping process on scale-free networks 
with modular structure remains less understood, in spite 
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of the facts that modularity plays an important role 
in shaping up scale- free networks [1^, and that taking 
into account the modular structure of scale-free networks 
leads to a better understanding of how the underlying 
systems work [SCj. 

In this paper, we study the classic trapping prob- 
lem on a class of hierarchical networks [3, Q, which 
is a random walk problem with a single immobile trap 
positioned at a given site, absorbing all walks visit- 
ing it. Here we focus on a particular case with the 
trap located at the node with the highest degree. The 
networks studied can capture simultaneously scale-free 
behavior and modular structure. Moreover, the net- 
works belong to a deterministic growing type of net- 
works, which have received much attention from the 
scientific communities and have proved to be a useful 
tool P, m, m, m m, [13, Hi, m H mi ■ The de- 
terministic nature of the hierarchical networks makes it 
possible to investigate analytically the trappi ng p rocess 
defined on them. By applying the formalism j62l . [63j of 
generating functions |64l | for random walks, we derive the 
rigorous solution to the mean first-passage time (MFPT) 
that characterizes the trapping process. The obtained 
exact result shows that the MFPT scales algebraically 
with the number of network nodes. We also compare the 
behavior of the trapping problem on the hierarchical net- 
works with those of other networks, and show that the 
hierarchical networks can be helpful for enhancing the 
efficiency of the trapping process. 



II. MODULAR SCALE-FREE NETWORKS 

Let us introduce the model for the hierarchical scale- 
free networks with a modular structure, which can be 
constructed in an iterative way 0, [1] . We denote by Hg 
the network model after g {g > 1) iterations (number of 
generations). Initially {g = 1), the network consists of 
a central node, called the hub (root) node, and M — 1 
peripheral (external) nodes with M > 3. All these ini- 
tial M nodes are fully connected to each other, forming 
a complete graph. At the second generation {g — 2), 
we generate M — 1 copies of Hi and connect the M — 1 
external nodes of each replica to the root of the original 
Hi. The hub of the original Hi and the {M ~ 1)^ periph- 
eral nodes in the replicas become the hub and peripheral 
nodes of H2, respectively. Suppose one has iJg-i, the 
next generation network Hg can be obtained from -ffg-i 
by adding M — 1 replicas of -ffg-i with their external 
nodes being linked to the hub of the original Hg^i unit. 
In Hg, its hub is the hub of the original ifg-i, and its 
external nodes are composed of all the peripheral nodes 
of the M — 1 copies of -ffg-i- Repeating indefinitely the 
replication and connection steps, we obtain the hierar- 
chical modular scale-free networks. Figure [T] illustrates 
the construction process of a network for the particular 
case of M = 5, showing the first three iterations. 

According to the network construction, one can see 




FIG. 1: The iterative construction process of a hierarchical 
network for the case of Af = 5. Notice that the diagonal nodes 
are also connected — links not visible. 

that Hg, the network of gth generation, is characterized 
by two parameters g and M , with the former being the 
number of generations, and the latter representing the 
replication factor. In Hg, the number of nodes, often 
called order of the network denoted as Ng , is Ng = . 
All these nodes can be classified into the following four 
sets [H, [6^: peripheral node set P, locally peripheral 
node set Pm (1 < < <?), set H only consisting of the 
hub node of Hg, and the local hub set IHI,„ (1 < m < g); 
see Fig. [5] The cardinalities, defined as the number of 
nodes in a set, of the four sets are 

\n = {M-ir, (1) 

|P„,| = (M - l)m]^p-{m+l)^ (2) 

|H| = 1, (3) 

and 

= (Af-l)M^'-("+i), (4) 

respectively. For Hg, all nodes belonging to the same 
set have identical connectivity (i.e., degree), which are 
known exactly. For example, the degree Kh{g) of the 
hub node is the largest; it has a value of 

Any node in P has the degree 

K^{g) = g + M-2. (6) 
Again, for instance, the degree of a node in P^ is 

Kp,m{g) = m + M-2, (7) 
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FIG. 2; (Color online) Classification of nodes in network 
for the case of M — 4. The filled circles, open circles, full 
square, and triangles represent peripheral nodes, locally pe- 
ripheral nodes, hub node, and locally hub nodes, respectively. 



FIG. 3: Labels of all nodes in H3 for the case of M 
corresponding to the g — 3 case in Fig. [T] 



and an arbitrary node in Mm has a degree of 



M - 1 

~ (8) 



M-2 

Thus, the sum of degrees for all nodes in Hg is 



features make the networks unique within the category 
of scale-free networks; it therefore is worthwhile to inves- 
tigate various dynamical processes running on them. In 
what follows we will study the trapping problem on this 
class of modular networks to uncover the influence of the 
particular topologies on the trapping process. 



Da 

9-1 9-1 
= Khig) + J2 KhMgmml + Kpig)\F\ + ^ A'p,,„(g)| 

m— 1 m— 1 

= {3M - 2){M - l)M<^~^ - 2{M ~ iy+\ 
and the mean degree averaged over all nodes is 



2D 

{k)g = ^ = 2{M-l) 



2 

3 2 

M \ M 



(10) 



which is approximately equal to 2(M — l)(3Af — 2) /M in 
the limit of infinite g. 

The hierarchical networks present some typical prop- 
erties of real systems in nature and society [65|. They 
are scale free with the degree distribution exponent 7 — 
1 -t- lnAf/ln(Af — 1). The average path length, defined 
as the shortest distance averaged over all pairs of nodes, 
scales logarithmically with the number of nodes. In the 
large network order limit, the average clustering coeffi- 
cient tends to a large constant dependent on Af . Thus, 
the whole family of networks exhibits small- world behav- 
ior [63] ■ In addition, the betweenness of nodes in the net- 
works follows the same power-law distribution Pb ~ 
irrespective of Af . In particular, the networks show an 
obvious modular structure. All these characteristics are 
not shared by other models. The peculiar topological 



III. FORMULATION OF THE TRAPPING 
PROBLEM 

- 7n I 

In this section we formulate the trapping problem on 
family of hierarchical scale- free networks i?g, which 
is actually a simple unbiased Markovian random walk of 
a particle in the presence of a trap or a perfect absorber 
located on a given node. To facilitate the description, 
we distinguish different nodes in Hg by assigning each 
of them a labeling in the following way. The hub node 
in Hg has label 1; the other Af — 1 peripheral nodes in 
Hi are labeled as 2, 3, and Af — 1, respectively. Assume 
that we have labeled nodes in ifg-i consecutively by 1, 2, 
and M^~-^; in the next generation g, we keep the labels of 
nodes in the original ffg-i unchanged and label only the 
nodes belonging to the Af — 1 copies of ifg-i by assigning 
to each node a different integer from M^~^ -I- 1 to . In 
this way, every node in Hg is labeled by a unique integer 
from 1 to Ng = MS; see Fig. E] 

For convenience, we continue to represent Hg by 
its adjacency matrix Ag of order iV^ x Ng, whose 

ttij 1 if 
— other- 



element Uij is defined as follows: 
i and j are neighboring nodes and a.y 
wise. Then the degree, di{g), of node i is given by 



d^{g) = E 



the diagonal degree matrix Zg of 



Hg is Zg 



diag(di(g), d2ig),...,d^{g),..., dNAg)), and 
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the normalized Laplacian matrix of Hg is provided by 
lig — Ig ~ Zg^ Ag, where Ig is the Ng x Ng identity 
matrix. 

Before proceeding further, let us introduce the so- 
called discrete-time random walk on Hg. At each time 
step, the particle jumps from its current location to any 
of its nearest neighbors with equal probability. Accord- 
ing to this rule, at time i, a particle located at a node 
i will hop to one of its di{g) neighbors, say u, with the 
transition probability aiu/di{g). Suppose that the par- 
ticle starts off from node i at t — Q, then the jumping 
probability Pij of going from i to i a, t time t is governed 
by the following master equation |22l |: 

We next focus the trapping problem on Hg with the 
trap fixed on the hub node, i.e., node 1, represented as 
It- The particular choice for the trap position allows to 
compute analytically the MFPT, which will be discussed 
in detail in the following section. Similar to the standard 
discrctc-time random walks, during the trapping process, 
in a single time step, the particle, starting from any node 
except the trap it, jumps to any of its nearest neighbors 
with the same probability. What we are concerned with is 
the expected time that the particle spends, starting from 
a source node before being trapped, which is in fact a 
random variable. Let Xi be the expected time, frequently 
called first-passage time (FPT) or trapping time, for a 
walker, starting from node i, to first arrive at the trap it- 
In order to determine Xi, we define F{Xi — t) to be the 
probability for the particle, starting from point i, to first 
hit the trap after t steps. Notice that since the Markov 
chain [68| representing such a random walk is ergodic, 
the particle will be eventually trapped independently of 
the origin, implying that X^t^o ^(^i — t) — 1 holds for 
all i. It is easily known that the set of these interesting 
quantities obeys the following recurrence relation: 

Fix, = t) = Y, -^MX^ =<-!), (12) 

where i ^ it- 

Let Fi{z) be the corresponding generating function of 
quantity F{Xi = t): 

oo 

= = i)z*, (13) 

t=o 

which encapsulates all the information contained in the 
discrete probability distribution F(Xi = t). For exam- 
ple, the expected value Xi is the first derivative of Fi{z) 
evaluated at z = 1. 

Let F(z) stand for the {Ng — l)-dimensional vec- 
tor [^2(2;), ^3(2;), . . . ,Fn^{z)]^, where the superscript T 



represents the transpose of the vector. According to 
Eqs. dn]) and we have 

F(z) = zWF(z), (14) 

where W is a matrix with order (Ng — 1) x (Ng — 1) with 
entry Wij — aij/di{g). Differentiating the two sides of 
Eq. pi]) with respect to z and doing some simple algebra 
operations, we have 

(I-zW)F'(z)- WF(z) = 0, (15) 

in which I is the identity matrix with order {Ng — 1) x 
{Ng — 1); is the (iV^ — l)-dimensional zero vector 
(0, 0, . . . , 0)^. Setting z = 1 in Eq. ^ leads to 

F'(l) = (I - W)-iWF(l) = (I - W)-ie, (16) 

where e = (l,l,...,l)^is the (iV^ — l)-dimensional unit 
vector. Actually, (I — W)~^ in Eq. is the funda- 
mental matrix of the Markov chain representing the un- 
biased random walk, and I — W is a submatrix of the 
normalized discrete Laplacian matrix Lg of Hg , which is 
obtained from Lg by removing from Lg the first row and 
column corresponding to the trap. 

From Eq. p6)) . the mean first-passage time, (T)g, 
which is the average of Xi over all initial nodes dis- 
tributed uniformly over nodes in Hg other than the trap, 
is given by 

" 1=2 i=2 j=2 

where kj is the corresponding {i,j) element of matrix 
(I — W)~^, which is the mean time that the particle 
spends at node i starting from node j 69] . 

Equation p7|) shows that the problem of calculating 
MFPT {T)g is reduced to finding the sum of all elements 
of matrix (I — W)^^. Since the order of (I — W) is 
{Ng — 1) x {Ng — 1), where Ng increases exponentially 
with g, for large g, inverting matrix (I — W) is pro- 
hibitively time and memory consuming, making it in- 
tractable to obtain (T)g through direct calculation from 
Eq. ([TT]); one can compute directly the MFPT only for 
the first several generations (see Fig. 0]). Hence, an al- 
ternative method of computing MFPT becomes neces- 
sary. In [t^I , to allow for a drastic reduction in computa- 
tional cost, a scheme was proposed mapping the original 
Markov process on another Markov process. Although 
the method can bring down the computational efforts, it 
is an approximate one. Fortunately, the special recursive 
construction of the hierarchical networks allows to cal- 
culate analytically MFPT to obtain an explicit solution 
for arbitrary generation g. In the next section, we will 
provide the detailed process for the derivation of MFPT 
using a method significantly different from that applied 

in gil,!!!!!!,!!!. 
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FIG. 4: (Color online) Mean first-passage time {T}g as a func- 
tion of the iteration g on a semilogarithmic scale for two cases 
of M — 3 and M = 4. The filled symbols are the data coming 
from genuine simulations of the trapping process; the empty 
squares and circles represent the numerical results obtained 
by direct calculation from Eq. (I17|l : while the empty pen- 
tagons and hexagons correspond to the exact values given by 

Eq. dsni)- 



IV. CLOSED-FORM SOLUTION TO MEAN 
FIRST-PASSAGE TIME 



Prior to deriving the general formula for MFPT, (T) 
for the trapping issue on Hg , we first define some related 
quantities. Let Pg{t) denote the probability that, at the 
generation g, the particle starting from any peripheral 
node in P first arrives at the hub after t jumps; and let 
Qg{t) represent the probability that, the walker originat- 
ing from the hub to first reach any node belonging to P 
after t steps. Then, the following fundamental relations 
can be established: 



'Jt.i 



M 



3-1 t-1 



1 



^^Q™«Fg(i-l-z) (18) 



the last term on the rhs describes the probability of the 
process in which the walker first makes a jump to a local 
hub node belonging to Pm, then it takes i time steps, 
starting from the local hub, to hit one of the nodes in P, 
and continues to jump more t — 1 — i steps to first reach 
the hub. 

Analogously, the two terms on the rhs of Eq. can 
be understood based on the following two processes. The 
first term explains the occurring probability of the pro- 
cess that the walker, starting from the hub, only needs 
one time step to reach a peripheral node in P. The sec- 
ond term represents the happening probability of such a 
process that the particle, originating from the hub, first 
makes one jump to a local peripheral node in Pm, then 
makes i jumps to the hub, and proceeds to the destina- 
tion (one of the nodes in P) , taking more t — 1 — i time 
steps. 

Equations (fT8|) and (fT9|l provide the two basic relations 
governing the trapping problem performing on Hg, from 
which almost all subsequent results are derived from. As 
shown in the preceding section, although we are con- 
cerned about only the fundamental quantity (namely, 
MFPT), the direct calculations are practically hard and 
intractable for large networks. Nevertheless, the particu- 
lar construction of the networks allows to overcome this 
difficulty in virtue of the powerful mathematical tech- 
nique of generating functions [64] , through which we can 
compute and determine the MFPT {F)g indirectly. 

First, we define two generating functions, Pg{x) and 
Qg{x), for the probability distribution of first-passage 
time described in Eqs. p8| and (fT9|) . which can be writ- 
ten as 



Pg{x) = X]P,(i)a 



M - 2 



XPg{x) 



X 

KM 



9-1 

^ Qm{x)Pg{x) 
m— 1 

(20) 



and 



and 



Qg{t) 



KM 



9-1 t-1 

-EE 

m— 1 -i— 1 



(M - If 
Kh{g) 



■P^{^Qg{t-l-i), 
(19) 



where 5t^i is the Kronecker delta function that is defined 
as follows: St^i — 1 \i t is equal to 1, and St.i — other- 
wise. Note that in Eqs. (fT5)) and the equivalence of 
nodes in the same set (e.g., P or P^,) was used. 

The three terms on the right-hand side (rhs) of Eq. (fT8|) 
can be elaborated as follows: the first term accounts for 
the probability that the walker takes only one time step 
to first reach the hub; the second term on the rhs explains 
the case that the particle gets first to one of its M — 2 
neighbors belonging to P in one time step, and then it 
takes more t — \ steps to first arrive at the target node; 



t=0 



(M - 1)9 
KM 



^MxYiM 

KM 



'Pmix). 

(21) 



After some algebraic operations, Eqs. pp]) and (PT|) can 
be recast, respectively, as 



PM 



Kpig) 



(M-2) 



9-1 



(22) 



6 



and 



and 



Qgix) 



Khjg) 1 

{AI - 1)9 X 



9-1 

^(M-ir-^F„(x) 



(23) 

Let denote the first-passage time for a walker start- 
ing from an arbitrary node in P to reach the hub for the 
first time, which is in fact the number of steps for the 
walker originating from any node in P to first visit the 
hub. Let stand for the FPT needed for a particle 
initially located at the hub to first hit any node in P. 
Then, according to the property of generating functions, 
the two quantities and are given separately by 



d ~ 
dx 



and 



T, 



H 



d ~ 
dx 



Qgix) 



(24) 



(25) 



x = l 



Differentiating, respectively, both sides of Eqs 
and (j23p with respect to x and setting a; = 1, we ob- 
tain the following two coupled relations: 



Af-2- 



9-1 



(26) 



and 



9 



Khjg) 1 

{M ~ 1)3 [M - 1)3 



9-1 



5](M-1)™T^. (27) 



From the above two coupled equations, it is not difficult 
to have 



T, 



9+1 



T. 



1 + T, 



H 



and 



(M-i)T«i-r/ = M-i + r/ 



(28) 



(29) 



Considering the initial conditions T2 — M + 1 and 
= {2M - 1)/(M - 1), we can solve the simultane- 
ous equations, i.e., Eqs. (|28p and (P^j) . to obtain 



= 3A/ - 8 



7M 



M2 



M 



M -1 



2M + 3 (30) 



H 



5M - 2 
M2 



M 



M -1 



(31) 



The obtained expressions for TF' and are very im- 



portant, using which we will determine MFPT {T)g. 



To 



facilitate the computation, we use Qg to represent the set 
of nodes in Hg and separate them into two subsets: one 
subset is f^g-i made up of nodes in the original Hg-i, and 
the other subset, denoted by Ctg, is the set of nodes of the 
AI—1 copies of Hg^i. Let Ti{g) denote the trapping time 
for a walker originating at node i on the gth generation 
network to first reach the trap node (hub). Obviously, 
for all (7 > 0, Ti{g) — 0. For 5 = 1, it is a trivial case, we 
have T2(l) = T^{1) = • • • = Tm{1) = M - 1. Then, by 
definition, the MFPT {T)g can be expressed as 



1 

3 1=2 



(32) 



where the sum term X]i='2 '^iid) be rewritten as 



Y^Ug) = Y.T,{g)+Y,Ug) 

= T.(g-l)-f ^r.(g), (33) 



1=2 



which is obvious from the particular construction of the 
hierarchical networks. Thus, we have 



(T) 



{T)g-, + -^yT,{g). (34) 



Hence, to obtain an exact solution for {T)g, all that is 
left is to evaluate the sum in Eq. ([34|) , with a goal to first 
find a recursive relation for {T)g. From Figs. [T] and [21 
the sum term on the rhs of Eq. (|34p can be evaluated as 
follows: 



Y Ug) = t;|P| + ^^(Tf + 1) + E (^-^ - 1)'"""' [(^™ - + ^™^™+i + ^^X] ■ (35) 

ieQg "1=1 

I 



Substituting previously obtained equations for the ex- pressions of related quantities in Eq. (|35p and combining 
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with Eq. (|33p. we can obtain the following recurrence 
relation for {T)g: 



{N,+,-l){T)g+,-M{Ng-l){T), 



2(M- 1) 
M2 



-MS 



(3M - 2) 



M 



M - 1 



7\/ 



(36) 



Using the initial condition (T)2 = M(Af + 2)/(Af + 1), 
Eq. ([36]) is solved inductively to obtain the rigorous ex- 
pression for the MFPT: 



{T)o 

M9-3(M- 1) 
M9 - 1 

f3 



(6^/=^ - IGM^ + 14M - 4) 



Af- 1 



(5Af3 - 10Af2 + 4Af) - 2.9 (Af ^ - Af ) 



(37) 



We have checked our analytic formula against numeri- 
cal values obtained according to the fundamental matrix 
provided by Eq. p7|) . For different parameters Af and 
g, the values obtained from Eq. ([57)1 completely agree 
with those numerical results on the basis of the direct 
calculation through Eq. (jlT)) : see Fig.[4l This agreement 
serves as an independent test of our theoretical formula. 
Moreover, we have performed genuine simulations of the 
random walk process on the hierarchical networks. The 
data from the true process are shown in Fig. [H each of 
which is obtained by averaging over 10 000 realizations. 
The results of the true simulations are in excellent agree- 
ment with our analytical ones given by Eq. (j37p . and 
thus provide an important further evidence in favor of 
our findings. 

We continue to show how to represent MFPT in terms 
of network order Ng with the aim to obtain the scaling 
between these two quantities. Recalling Ng = , we 
have g = log^^ Ng . Hence, Eq. ([57)1 can be rewritten as 



M - 1 



N„ 



AP 



Ng-1 



(6Af^ - 16Af2 + 14A/ - 4) {Ng 



those scalings for other media. For instance, on regular 
lattices with large order N, the asymptotical behavior of 
MFPT (T) is (T) - N^, (T) - iVlniV, and (T) - N for 
dimensions d — 1, d — 2, and d — 3, respectively |36|. 
Again for example, on planar Sierpinski gasket [37| and 
Sierpinski tower [38l| in three Euclidean dimensions, and 
the T fractal HjTthe MFPT (T) scales superlinearly 
with network order; i.e., it grows as a power-law func- 
tion of network order with the exponents being 1.464, 
1.293, and 1.631, respectively. Fi nally , for the pseud- 
ofractal web [4l|, the Koch network [i^] , and the Apollo- 
nian network [43|, they are all scale free, their MFPT 
scales linearly or sublinearly with network order, fol- 
lowing separately the asymptotical behaviors (T) ^ N, 

(T) - Arln2/ln3^ ^^^^ (j^^ ^ ^2-ln2/l„5^ rpj^^g^ 

pared with the aforementioned regular networks, frac- 
tals, even scale-free networks, the addressed hierarchical 
networks exhibit more efficient configuration for random 
walks with a single trap fixed at the node with highest 
degree. 

The root of the high efficiency of the trapping prob- 
lem on the hierarchical scale-free networks lies in their 
architecture. In this network family, there are many 
small densely interconnected clusters, which combine to 
form larger but less compact groups connected by nodes 
with high degrees (i.e., local hub nodes). The relatively 
large groups are further joined to shape even larger and 
even less densely interlinked modules. These modules or 
groups are combined again at a "large" node forming a 
fine hierarchical structure that is responsible for the fast 
diffusion phenomenon, which can be understood from the 
following heuristic argument. When a walker starts off 
from some node, it will either hit the hub directly or 
first get to local hub nodes. These local hubs, although 
not connected to the trap node, play the role of bridges 
linking different modules together at the local peripheral 
nodes, through which the walker may easily find the way 
to the trap. Thus, the walker can visit the trap in a short 
time, disregarding its starting points. 



-(5A//3 - lOM^ + AM) - 2 (Af2 _ m) log^^ Ng 



(38) 



V. CONCLUSIONS 



Thus, for networks with large order, i.e., Ng 



CXD, 



(Ng^W = {Ng] 



l-ln(M-l)/lnM 



(39) 



where the exponent d{M) is lower than 1. Clearly, 9{M) 
is a decreasing function of M: when Af grows from 3 to 
infinite, 6'(Af) descends from 1 — In 2/ In 3 and approaches 
to zero, which means that the efficiency of the trapping 
process depends on Af . The larger the value of Af , the 
more efficient the trapping process. Equation ([55]) also 
implies that in the infinite network order Ng limit, the 
MFPT grows algebraically with increasing order of the 
networks. 

The above obtained scaling of MFPT with order of 
the hierarchical scale- free networks is quite different from 



In conclusion, we have investigated the classic trapping 
problem on a class of hierarchical networks that can bring 
under a single roof the scale- free and modular topologies, 
which are two striking structural properties observed in 
various biological and social networks. Thus, the hierar- 
chical networks can mimic some real-world natural and 
social systems to some extent (to what extent it does is 
still an open question). Using the method of generat- 
ing functions, we derived the recursion relations govern- 
ing the evolution of the MFPT for random walks on the 
networks, with the only trap located at the hub node. 
These recursive relations are obtained from the special 
construction of the networks, from which we determined 
explicitly the solution for the MFPT, which shows that 
the MFPT {T)g varies algebraically with network order 
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Ng as {T)g ~ (iVg)^(^^) wltL the exponent 0{M) much 
less than 1 that decreases from 1 — In 2/ In 3 to zero when 
M increases from 3 to infinite. Thus, in the fuU range 
of M, the efficiency of the trapping process on the hi- 
erarchical networks is high. We have also compared the 
result with those previously obtained for other media, 
and found that in marked contrast to other graphs, the 
hierarchical networks have more efficient structure that 
tends to speed up the diffusion process. Finally, it de- 
serves to be mentioned that although the hierarchical 
networks are efficient for the trapping problem with the 
trap fixed on the hub, they might lose this characteris- 
tic when the trap is positioned at a randomly selected 
node, due to the somewhat tree-like macro-structure of 
this kind of networks [1^ . 
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